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We present a study of the X(3940) state in the process e + e~ — > J/tp D* D. The X(3940) mass and 
width are measured to be (3942 ± 7 6 ± 6) MeV/c 2 and T = (37 1 ^ ± 8) MeV. In the process e+e" -> 
J/tp D* + D*~ we have observed another charmonium-like state, which we denote as X(4160), in 
the spectrum of invariant masses of D* + D*~ combinations. The X(4160) parameters are M = 
(4156 1 20 ± 15 ) MeV/c 2 and Y = (139 ± ± 21) MeV. The analysis is based on a data sample with 
an integrated luminosity of 693 fb -1 recorded near the T(4S) resonance with the Belle detector at 
the KEKB e + e~ asymmetric-energy collider. 

PACS numbers: 13.66.Bc,12.38.Bx,14.40.Gx 



Double charmonium production in e + e~ annihilation, first observed by Belle in 2002 [1], can be used to search for 
new charmonium states, recoilingagainst some known and easily reconstructed charmonium. The study of various 
double charmonium final states [2|, 13J demonstrated that there is no significant suppression of the production of 
radially excited states: the cross-sections for J/tprj c , tp(2S)r] c , J/tpr] c (2S) and tp(2S)r/ c (2S) are very close to each 
other. These studies also show that scalar and pseudoscalar charmonia are produced copiously recoiling against J/tp 
or ■0(25'). A new charmonium-like state, X(3940), has been already observed in the spectrum recoiling against J/tp, 
and reconstructed in the D*D [4j final state [5J. On the other hand, there has recently been a number of reports 
on observation of new charmonium or charmonium-like states above DD threshold [6]. Their properties are quite 
different from those expected from the quark model. These experimental results have renewed theoretical interest in 
spectroscopy, decay and production of charmonia |7|. 

In this Letter we present a new study of the X(3940) resonance and report on the observation of a new charmonium- 
like state in the process e + e~ — > J/i/j D^*'D^ and the measurement of its parameters. The integrated luminosity used 
for this analysis is 693 fb -1 collected with the Belle detector Q near the T(45) resonance at the KEKB asymmetric- 
energy e + e~ collider 0. 

This study is performed using the selection procedure similar to that described in Ref. 0, [B| . All charged tracks 
are required to be consistent with originating from the interaction point. Charged kaon candidates are required to 
be positively identified, while no identification requirements are applied for pion candidates as the pion multiplicity 
is much higher than those of other hadrons. Kg candidates are reconstructed by combining 7r + 7r~ pairs with an 
invariant mass within 10MeV/c 2 of the nominal Kg mass. We require the distance between the pion tracks at 
the Kg vertex to be less than 1 cm, the transverse flight distance from the interaction point to be greater than 
1mm and the angle between the Kg momentum direction and decay path to be smaller than 0.1 rad. Photons are 
reconstructed in the electromagnetic calorimeter as showers with an energy more than 20 MeV that are not associated 
with charged tracks. Photons of energy more than 50 MeV are combined to form tt° candidates. If the mass of 77 
pairs lies within 15MeV/c 2 of the nominal 7r° mass, such pairs are fitted with a 7r° mass constraint and considered 
as 7T° candidates. J/tp candidates are reconstructed via the J/tp^l + £~ (I — e, /i) decay channel. Two positively 
identified lepton candidates are required to form a common vertex that is less than 1 mm (« 6 a) from the interaction 
point in the plane perpendicular to the beam axis. A partial correction for final state radiation and bremsstrahhmg 
energy loss is performed by including the four-momentum of every photon detected within a 50 mrad cone around 
the electron direction in the e + e~ invariant mass calculation. The J /tp signal region is defined by the mass window 
\Mi+£- — Mj/^ J < 30 MeV/c 2 (w2.5 a). J/tp candidates are subjected to a mass- vertex fit to improve their momentum 
resolution. QED processes are suppressed by requiring the total charged multiplicity in the event to be more than 4. 
J/tp mesons from BB events are removed by requiring a center-of-mass (CM) momentum V*ji^ > 2.0GeV/c. 
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We reconstruct D° mesons using five decay modes: K 7r+, K K + , K it n + ir + , KgTr + n and K 7r + 7r°. Candi- 
date D + mesons are reconstructed using K~ir + ir + , K~K + -k + and K%-k + decay modes. A ±15MeV/c 2 mass window 
is used for all modes except D° — > K~n + ir° (±20 MeV/c 2 ) (~ 2.5 a in each case). To improve their momentum 
resolution, D candidates are refitted to the nominal D° or D + masses. To study the contribution of combinatorial 
background under the D peak, we use D sidebands selected from a mass window four times as large. For the study 
of the process e+e~ — * J/ipD*D^ > we use only the cleanest D* + — s-D°7r+ channel. D* + candidates from the signal 
window, selected in the interval ±3MeV/c 2 of the nominal D* + mass (« 2.5 c), are refitted to the nominal D* + 
mass. The D*+ sideband region is defined by 2.016 GcV/c 2 < M(D°ir+) < 2.028 GeV/c 2 . Only one J/ipD or one 
J /4>D* + combination per event is accepted; the combination with the best sum of \ 2 01 the mass fits for J/ip and 
candidates is selected. In the sidebands a single candidate per event is selected as well. The sideband is 
divided into windows of the same width as the signal one, and the candidate with the smallest difference in mass from 
the center of its window is chosen. 

The method for reconstructing the processes was described in [5j. In addition to the fully 

reconstructed J/tp, only one of the £)W's is fully reconstructed (referred to below as D^l- D xec — D° or D + , 
D* CC = D* + ), and the other unreconstructed (referred to as associated D = fisLc) in the event is observed as a 
peak in the spectra of masses recoiling against the reconstructed combination J/ ipD^ . The recoil mass against the 
particle or combination of particles is defined as 

M Iecoil (X) = \/{E C M -E* x f-p* x 2 , (1) 

where E x and p x are the CM energy and momentum of the (combination of) particle(s). The Af rcco ii( J/ipDicl) peak 
around the nominal mass of I?isLc with a typical resolution ~30MeV/c 2 is used to identify the studied process. As 
the resolution is smaller than Mr>* — Mjj, the method allows the contributions from the processes e+e" — > J/tpDD, 
J /tj}D*D and D*D* to be disentangled. The M r0C oii( J/ipD rcc ) and M Teco n(J/ipD* ec ) spectra in the data are shown in 
Fig. [T] as points with error bars for the signal Dree windows; histograms show the scaled Dree sideband distributions. 
The signals for the processes e + e~ — > J/ipDD, D*D and D*D* are evident in Fig. Q]a) at the D and D* nominal 
masses and at a mass ~2.2GcV/c 2 , respectively. The latter peak is shifted and widened due to two missing pions 
(or photons) from D* decays. Another excess at ~2.45GeV/c 2 can be explained by the process e + e~ — * J/ipDD**. 
The processes e + e~ — > J /ipD*D and D*D* are also clearly seen in Fig. Q]b) as distinct peaks around the D and D* 
nominal masses. We use -D re c sidebands to describe the combinatorial background contribution through simultaneous 
likelihood fits to the D^l oc signal and sideband spectra. The signal shapes are fixed from the Monte Carlo (MC) 
simulation. The background distribution is parameterized by a second-order polynomial function (linear function in 
case of D* ec ). Only the region below 2.35GeV/c 2 is used because of a possible contribution from e + e~ —> J/ipDD**. 
The signal yields (including the tail due to initial state radiation [ISR]) and statistical significances are listed in 
Table D 

TABLE I: Summary of the results of the fits to the M leco n(J/ipD lac ) and M reC oii( J/tl>D* ec ) spectra. 





J/^D ICC 


J/i>D*„ 


JC 




N 


Ma 


N 


Ma 


e+e" ^J/i>DD 


162 ± 25 


7.6 






e+e" ^J/ipD*D 


159 ±28 


6.5 


19.0 +11 


5.8 


e+e" ^J/ipD*D* 


173 ± 32 


5.6 


47.2 1?;! 


8.4 



We perform a study of these observed processes and search for new charmonium states X c5 that can be produced via 
e+e" -^J/ipX c5 followed by the decay X cS ^ . Tagging the process e+e" — * J/tjj £)(*)£)(*) by the requirement 

\M reca {\{J / tpD^c) — M D (,) | < 70 MeV/c 2 we thus divide each of selected J/ipD or J/ipD* + combinations into two non- 
overlapping samples, each comprising ~ 50% of the signal events. The ISR tail causes an inefficiency for the tagging 
requirement as well as cross talk between different final states: the contribution of the process e+e" — > J/ipDD (D*D) 
to the sample tagged as e+e" — > J/ipD*D (D*D*) is ~ 10%, the reverse cross talk is only ~ 1.5% and neglected. In 

our study we constrain M Ieco n(J/ipDTec) to the D^ s l oc nominal mass. This improves the resolution on M Teco n(J / , 
which corresponds to the invariant mass of the produced D^*' meson pair, by a factor of 3 — 10 with respect to the 
unconstrained value (~ 30MeV/c 2 ). The M(Z)W5>)) reso lution varies from 2MeV/c 2 at threshold to 8MeV/c 2 
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FIG. 1: The distributions of masses recoiling against the reconstructed a) J/ipD and b) J/ipD* combinations in the data. 
The histograms show the scaled D^*' sideband distributions. The solid curves are results of the fit, the dashed curves are the 
background functions. 



at M(£>WZ)W) _ 5.0GeV/c 2 for all the processes except e + e'^J/^D*D with Aec Assoc- In the latter case the 
resolution is worse because of the Acc from the D* decay (-10MeV/c 2 at M(£>*A)~3.94GeV/c 2 ). 
In the data the spectra of M(flWl)W) are shown in Figs.[2]a), b), c), d) for A C c Assoc, A 



D* D* 

rcc assoc 



cases, respectively. Points with error bars correspond to the D 



(*) 



D* D* D 

^ assoc , -^rcc-^assocj 



signal windows while hatched histograms 
show the scaled Acc sideband distributions. Excesses from the signal Acc window over the sideband distributions 

-> D*D*) estimated using the 
processes are generated with 



D*D and D*_D 



are seen around the threshold in all figures. The reflections (DD 
MC simulation are shown with open histograms. In the MC the e + 
MpW5(*') spectra tuned to the data. 

We perform simultaneous likelihood fits to Aec signal and sideband distributions to fix the combinatorial back- 
ground shapes. The accuracy of description of combinatorial backgrounds by Acc sidebands is validated with the MC 
simulation and with the data, where the M( A*) A*) ) spectra in the different sideband intervals are found to be in good 
agreement with each other. The combinatorial backgrounds are parameterized by the function A^/M — A/ t h r ■ e~ B ' M , 
where A and B are free parameters, except for the case A™ Assoc, where this shape is found to describe poorly the 
behavior of the background. In the latter case we parameterize the background by a relativistic Breit-Wigner function 
with a free mass, width and amplitude. The signal functions are a sum of a relativistic s-wave Breit-Wigner function 
and a threshold function (\/M — Mthr) to account for possible non-resonant production. The signal functions are 
convolved with the resolution functions and multiplied by the efficiency function obtained from the MC simulation. 
The reflections are taken into account in the fit. 

The fitted parameters of the Breit-Wigner functions and significances of the resonance contributions are listed 
in Table ITT1 We assess the significance of each signal using — 21n(£ /'Cmax)j where £ max is the maximum likelihood 
returned by the fit, and Cq is the likelihood with the amplitude of the Breit-Wigner function set to zero. This quantity 
should be distributed as x 2 ( n dof = 3) in the absence of signal, as three signal parameters are free in the fit for £ max - 
The non-resonant contributions are consistent with zero within la in all fits, except for the case A*cc^assoc (1.6 cr 
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FIG. 2: The M (-Dree -DasLc) spectra for events tagged and constrained as a) e + e —*J/ipDD, b), c) e + e — » J /ipD* D and d) 
e + e~ — > J /ipD* D* in the data. 

from zero). The fit results are shown in Fig. [2] as the solid curves; the dashed curves are the background functions. 
The insets in Fig. [2]a) and b) show the background subtracted spectra with the signal functions superimposed. 

A fit to M(DD) distribution finds a broad resonance near the threshold, which is tentatively denoted as X(3880), 
with a statistical significance of 3.8 er. However, the fit is not stable under variation of background parameterization as 
well as variation of the bin width. The fit with two resonances better describes the spectrum and is more stable, but 
the significance of the second resonance is lower than 3 a. We conclude that the observed threshold enhancement is not 
consistent with non-resonant e + e _ — * J/ipDD production, but with the present statistics the resonant structure in this 
process cannot be reliably determined. The significance of the X(3940) signal found by the fit to the M(Z? rcc Z)* ssoc ) 
spectrum is 6.0 c. The fitted width of X(3940) is slightly higher than that obtained in our previous analysis The 
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TABLE II: Summary of the signal yields, masses [MeV/c 2 ], widths [MeV] and significances for e + e — > J/ip(D^*' D^-*') 



State 


-AT events 


M 


r 


Ma 


X(3880)(D rcc O assoc ) 
X(3940)(L» rcc D: sS oc) 
X(3940)(Z) r * CCJ D asS oc) 
X(4160)(D* CC DL SOC ) 


63 

52 tli 

5.212-* 

no o + 12.3 
23.8 _ g0 


3878 ± 48 
3942 1 1 
3934 1?? 
4156 l 2 2 f, 


047 +316 
oq: ' - 143 
07 + 26 
°' - 15 
c 7 + 62 
°< -34 

1391 1 " 


3.8 
6.0 
2.8 
5.5 



TABLE III: Summary of the systematic errors in the masses (M in MeV/c 2 ), widths (T in MeV) and production cross sections 
(a in %) for X (3940) [X(4160)] resonances. 



Source 


X(3940) 


A(4160) 


M 


r 


a 


M 


r 


a 


Fitting procedure 


±4 


±6 


±5 


±12 


±18 


±2 


Selection 


±4 


±5 


±4 


±8 


±11 


±5 


Momentum scale 


±3 






±3 






Angular distributions 






±12 






±16 


Reconstruction 






±6 






±8 


Identification 






±4 






±4 








±3 






±4 


Total 


±6 


±8 


±16 


±15 


±21 


±20 



mass of the state is in good agreement with the reported mass, and the signal yield scales with respect to the previous 
result in proportion to the luminosity. Separate fits to the D® cc and D+ c samples yield 42 1 ^ and 8 1 \ signal events, 
respectively, in good agreement with the MC expectations (40 and 12) normalized to the integrated yield assuming 
equal branching fractions of A" (3940) decays into charged and neutral DD pairs. The X(3940) signal is also seen 
in the M (D* CC D assoc ) spectrum with a significance of 2.8 a, with parameters in good agreement with those from the 
M (_D rcc Z)* ssoc ) fit. As this sample is a small subsample of the D roc D* ssoc case, we use the latter fit for only as a cross 
check. The M(D*D*) spectrum demonstrates a clear broad enhancement around the threshold, which we denote as 
AT (4160). The AT (4160) signal is seen above the small combinatorial background and the X(3940) reflection with a 
statistical significance of 5.5 a. 

The Born cross sections for the processes e+e~ -> J/^AT(3940) [AT(4160)] multiplied by B Dl , m , = B(X -> D ( *W*) 
are calculated from the fitted X(3940) and A" (4160) yields with the procedure used in the previous analysis [2]. Taking 
into account the reconstruction efficiencies obtained from the MC simulation, the calculated Born cross-sections are: 

a(e+e- -» J/ipX(3940))B D ,^ = (13.9±|;t ± 2 - 2 ) * 
cr(e+e- -> Jli>X{Aim))B Dt p, = (24.7± 1 |j ± 5.0) fb. (2) 

The systematic errors of the parameters and production cross sections for AT(3940) and AT (4160) resonances are 
summarized in Table [TTTT To estimate the fitting systematics we study the difference in AT(3940) [AT (4160)] parameters 
returned by the fit to the Fig.[5]b) and d) distributions under variation of the signal and background parameterizations, 
the fit ranges and the histogram bins as well as the resolution functions. We also vary the definitions of the signal 
and sideband regions to check the stability of the resonance parameters. Another uncertainty in the determination 
of the masses is due to possible momentum scale bias. This was estimated in the previous paper ,5] to be smaller 
than 3MeV/c 2 . The systematic error for the cross section calculation is dominated by the uncertainty in the J/tp 
production and polarization angular distributions. In the MC both angular distributions are assumed to be flat and 
extreme cases (1 ± cos 2 9 and sin 2 9) are considered to estimate the systematic uncertainty in this assumption. In the 
case of the AT (4160) another source of the systematic uncertainty is the D* + polarization, which is also taken into 
account by varying the D* + helicity angle distribution. Other contributions come from the uncertainty in the track 
and 7T° reconstruction efficiencies; lepton and kaon identification and in the absolute B{D^). 

In summary, we have observed the processes e + e~ — > J/ijjDD (D*D, D*D*) and found significant enhancements 
in M(D^ £)W) spectra around thresholds in all these processes. A broad enhancement in M{DD) is not consistent 
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with non-resonant e + e~ —> J/t/jDD production, however the present sample is not large enough to allow the resonant 
structure in this process to be determined. We have confirmed our observation of the charmonium state, X(3940)— > 
DD*, produced in the process e + e~ J / ip X (39A0) with a significance of 5.7 a including systematics. The A(3940) 
mass and width are (3942 lg ± 6) MeV/c 2 and T = (37 ± 8) MeV. These measurements are consistent with our 
published results and supersede them. In this study we have found that the inclusive peak in the M reco n(J/ip) spectrum 
may consist of several states, thus our previous measurement of A (3940) branching fractions may be not reliable [f|. 
We report observation of a new charmonium-like state the A (4160) in the processes e + e~ — * J/ip X(4160) decaying 
into D*D* with a significance of 5.1ct, including the systematic uncertainty of the fit. The A (4160) parameters are 
M = (4156 1 20 ± 15)MeV/c 2 and T = (1391 n e { ± 21) MeV. 
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